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Abstract 

Background: Leishmaniasis is a group of diseases with a variety of clinical manifestations. The form of the disease 
is highly dependent on the infective Leishmania species and the immunological status of the host. The infectivity of 
the parasite strain also plays an important role in the progression of the infection. The aim of this work is to 
understand the influence of the natural infectivity of Leishmania strains in the outcome of visceral leishmaniasis. 

Methods: In this study we have characterized four strains of L infantum in terms of molecular typing, in vitro 
cultivation and differentiation. Two strains were isolated from HIV+ patients with visceral leishmaniasis (Bibiano and 
E390M), one strain was isolated from a cutaneous lesion in an immunocompetent patient (HL) and another internal 
reference strain causative of visceral leishmaniasis (ST) also from an immunocompetent patient was used for 
comparison. For this objective, we have compared their virulence by in vitro and in vivo infectivity in a murine 
model of visceral leishmaniasis. 

Results: Molecular typing unraveled a new k26 sequence attributed to MON-284 zymodeme and allowed the 
generation of a molecular signature for the identification of each strain. In vitro cultivation enabled the production 
of promastigotes with comparable growth curves and metacyclogenesis development. The HL strain was the most 
infective, showing the highest parasite loads in vitro that were corroborated with the in vivo assays, 6 weeks post- 
infection in BALB/c mice. The two strains isolated from HIV+ patients, both belonging to two different zymodemes, 
revealed different kinetics of infection. 

Conclusion: Differences in in vitro and in vivo infectivity found in the murine model were then attributed to 
intrinsic characteristics of each strain. This work is supported by other studies that present the parasite's inherent 
features as factors for the multiplicity of clinical manifestations and severity of leishmaniasis. 

Keywords: Leishmania infantum, Clinical isolates, Visceral leishmaniasis, Molecular typing, Metacyclogenesis, 
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Background 

Parasites from the Leishmania genus are trypanosomatid 
protozoans responsible for a group of diseases with a 
broad range of clinical manifestations collectively known 
as leishmaniasis (reviewed in [1-3]). The emergence of 
leishmaniasis as an opportunistic infection in HIV+ pa- 
tients in areas where both pathogens are endemic [4] 
has generated new interest in leishmaniasis. 

It is well known that species such as L. major and L. 
mexicana are usually exclusively dermotropic, while L. 
infantum and L. donovani are responsible for both cutane- 
ous and visceral leishmaniasis [5]. Apart from a general 
species-specific organ tropism of Leishmania, intraspecies 
intrinsic characteristics are also a relevant factor to con- 
sider. According to Maia et al. [6], dermotropic and 
viscerotropic L. infantum strains modulate the sand fly 
biting time on the host leading to the delivery, respect- 
ively, of a high or low dose of metacyclic promastigotes 
into the skin which will impact on the parasite tropism 
and manifestation of the disease. Even strains belonging 
to the same zymodeme have been associated to differen- 
tial infectivity [7] . 

In experimental infections, however, another parasite- 
related feature is of major importance. In vitro cultiva- 
tion of Leishmania is a subject open to wide variation 
between laboratories, making the comparison of similar 
experiments ambiguous. Depending on the culture 
medium (Santarem, N. and Cunha, J., submitted results 
and [8]), the duration of the culture [9] and the number 
of axenic passages performed [9], the promastigotes gen- 
erated will be differentially enriched in metacyclic forms 
[9], which will condition the success of the infection. 
Nonetheless, the genetics and the immune status of the 
host play a similarly important role in the tropism and 
severity of the disease [10]. In the murine models, L. 
major was only found in the infection site of the resist- 
ant C57BL/6 mice after subcutaneous injection, whereas 
the same experimental protocol followed in the suscep- 
tible BALB/c strain allowed visceralization [11]. Also, 
high and low infective strains maintained their profile 
(visceralizing or regulatory, respectively) in BALB/c and 
C.B-17 SCID mice, although with higher parasite loads 
in the T and B cell-dysfunctional SCID animals [12], 

The analysis of HIV/Le/s/zra<2m<2-coinfected human pa- 
tients brought important insights into the role of the im- 
mune system on the severity of the disease. On the one 
hand the visceralization of dermotropic strains is fre- 
quently observed in I Leishmania- co'mtections [13], 
as well as the regular presence of amastigotes in uncom- 
mon locations such as the lungs or the intestine [14]. 
On the other, the appearance of unique Leishmania 
zymodemes in HIV+ patients has been reported, which 
may be indicative of circulating strains normally associ- 
ated with asymptomatic disease in immunocompetent 



patients [13,15]. Some studies have shown that strains 
originating from HIV+ patients have low infectivity, 
which explains its appearance only in immunocomprom- 
ised individuals [7,16]. On the contrary, three distinct in- 
fective profiles were attributed to strains responsible for 
CL or VL (from immunocompetent or HIV+ patients) and 
no correlation was made according to the origin of the 
isolate [17]. 

In this study, we have focused on four different L. 
infantum strains isolated from patients with CL, VL and 
HIV '/ 'Leishmania coinfections. We characterized these 
strains according to molecular, biological and infectivity 
characteristics. We standardized the in vitro culture to 
avoid any biased infectivity that was evaluated with 
macrophage and mouse models. We have studied the 
distribution of the strains in acute and chronic infection 
by qPCR assessing the parasite load in spleen, liver, 
bone marrow, blood and lymph nodes and correlated 
differences in infectivity with major findings on the mo- 
lecular typing. 

Methods 

Parasites 

Four L. infantum strains isolated from patients in the 
Mediterranean basin and Portugal were used in this study. 
MHOM/MA/67/ITMAP-263 (ST) is a cloned line derived 
from a patient with visceral leishmaniasis [9,18] that was 
used as an internal and comparative control in all the ex- 
periments performed. HL strain (MHOM/PT/2009/LLM- 
1708) was isolated from an immunocompetent patient 
with cutaneous leishmaniasis. Briefly, a skin biopsy was 
dissociated in a cell strainer to isolate the cells and was 
then transferred into culture in RPMI. E390M (MHOM/ 
ES/99/LLM-855) and Bibiano (MHOM/ES/01/LLM-1083) 
were isolated from bone marrow aspirates of HW/Leish- 
mania-coinfected patients, this second one being respon- 
sible for recurrent relapses of leishmaniasis. The bone 
marrow samples were cultivated in NNN medium at 26- 
27°C, until the expansion of the promastigotes. ST, HL 
and Bibiano strains have been characterized by multilocus 
enzyme electrophoresis (MLEE) as MON-1 zymodeme, 
while E390M is MON-284 (electrophoretic mobilities for 
malate dehydrogenase (MDH) and glucose-phosphate 
isomerase (GPI) were determined to be of 104 and 105, 
respectively, in relation to MON-1 zymodeme [13]). 

Molecular typing 

L. infantum isolates were subjected to molecular typing 
by targeting four different regions of the Leishmania 
genome. Leishmania DNA was extracted by phenol/ 
chloroform as described below in more detail and sam- 
ples were adjusted to a final concentration of 10 ng/|iL 
after measuring DNA content with a Nanodrop ND- 
1000 spectrophotometer (Thermo Scientific). A volume 
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of 5 \iL of each sample was used in further PCRs. PCR 
products were run on 2% agarose gels stained with eth- 
idium bromide and visualized under UV light. Then they 
were excised from agarose gels and purified using the 
QIAquick Gel Extraction Kit (QIAGEN). First, the spe- 
cies status of the isolates was confirmed by DNA se- 
quencing of the heat-shock protein 70 (hsp70) gene [19]. 
Further subtyping was performed by sequence analysis 
of the ribosomal internal transcribed spacer 1 (ITS1) 
and 2 (ITS2) [20] and the hydrophilic acylated surface 
protein B (haspB) or I<26 gene [21]. The Big-Dye Ter- 
minator Cycle Sequencing Ready Reaction Kit V3.1 and 
the automated ABI PRISM 377 DNA sequencer (Applied 
Biosystems) were used for direct sequencing of the I<26, 
ITS1 and ITS2 PCR products that was performed with 
the corresponding forward and reverse primers; internal 
primers for sequencing were also used for the hsp70 
PCR product, as described by Fraga et al. [19]. The 
obtained sequences were analyzed and edited using the 
software BioEdit Sequence Alignment Editor, version 
7.0.9.0 (Ibis Biosciences) [22]. ClustalW multiple align- 
ment algorithm tool and manual adjustment were used 
for comparison of the resulting sequences with the re- 
spective published sequences. The hsp70 sequences were 
compared with those of the different Leishmania species 
generated by Fraga et al. [19]. ITS types were assigned 
to each isolate according to the sequence polymorphism 
of the 12 microsatellite regions included in ITS1 (four 
sites) and ITS2 (eight sites), as described by Kuhls et al. 
[20]. I<26 genotypes were assigned according to the size 
and sequence of the PCR product, following the criteria 
previously described by Haralambous et al. [21]. 

For the generation of unique patterns that could be 
used for strain identification, we amplified a region of 
the kinetoplast DNA minicircles and evaluated the re- 
striction profile after Haelll (Roche Applied Science) 
endonuclease digestion [23] using DNA extracted from 
axenic promastigotes and from experimentally infected 
murine tissues. 

Culture media 

Novy-MacNeal-Nicolle medium (NNN) was prepared with 
a semi-solid phase made of 1.4% agar (Sigma-Aldrich), 
0.6% NaCl (Merck), 31% defibrinated rabbit blood, 625 
units/mL penicillin, 625 units/mL streptomycin and RPMI 
1640 medium supplemented with 10% fetal bovine serum 
(FBS), 2 mM L-glutamine, 100 U/mL penicillin, 100 U/mL 
streptomycin and 20 mM HEPES buffer (all from Lonza) 
as liquid phase. 

Growth curves and viability 

Parasites were first passed in mice to control their viru- 
lence and frozen in vials for future use until 10 in vitro 
passages [9]. Promastigotes were cultivated at 26°C with 



an initial inoculum of 10 6 parasites/mL in NNN from a 
synchronized culture in the same media and followed 
for 6 days. In each day, parasites were counted in a 
hemocytometer and stained with Annexin V and 7- 
amino-actinomycin D (7-AAD) for viability analysis as 
described in [9]. 10 000 gated events were analyzed in a 
FACSCanto II (BD Biosciences) and the percentage of 
Annexin77AAD" cells determined with Flowjo software 
(TreeStar). 

Cell cycle 

In each day of culture, promastigotes were recovered and 
washed in PBS/FBS 2%. 2xl0 6 parasites were resuspended 
in 1 mL of PBS/FBS 2% and 3 mL of ice-cold absolute 
ethanol (Panreac) was carefully added while vortexing. 
After fixation for 1 hour at 4°C, the parasites were washed 
in PBS and resuspended in 1 mL of propidium iodide (PI) 
staining solution consisting of citrate buffer 3.8 mM in 
PBS, 50 (ig/mL PI (Sigma-Aldrich) and 0.5 \igl\iL RNAse 
A (Sigma-Aldrich). Following an incubation of 30 minutes 
at 4°C, 20 000 single live cells were acquired in a 
FACSCanto II and analyzed with the Flowjo s cell cycle 
built-in tool. 

Quantification of metacyclogenesis-dependent gene 
transcription 

10 7 promastigotes from day 1 to 6 of culture were 
resuspended in TRIzol reagent (Invitrogen) and frozen 
at -80°C. Total RNA was extracted using chloroform 
and isopropanol according to the manufacturers in- 
structions and solubilized in 10 \iL of nuclease-free 
water. RNA of high quality was obtained (RQI between 
9.0 and 10.0) as assessed using RNA StdSens Chips of 
the Experion automated electrophoresis system (Bio- 
Rad). RNA concentration was determined using a 
Nanodrop ND-1000 spectrophotometer. Samples were 
stored at -80°C until cDNA was synthesised. Reverse 
transcription was performed with iScript cDNA synthesis 
kit (Bio-Rad) according to the manufacturer s instructions 
over 500 ng of total RNA. Meta-1, Small Hydrophilic 
Endoplasmic Reticulum-associated Protein (SHERP) and 
histone H4 transcription was quantitatively analyzed after 
normalization with rRNA45 transcription by qPCR using 
the iQ SYBR Green Supermix according to the manufac- 
turers instructions in a My iCycler iQ5 (Bio-Rad). 4 |iL of 
cDNA (diluted 25 x) was used as template that was run in 
duplicate with 500 nM (Meta-1 and histone H4) or 250 
nM (SHERP and rRNA45) of the following primers (from 
Stabvida): Meta-1 [GenBank: NC_009401] forward: 5'- 
GGGCAGCGACGACCTGAT-3 ' and reverse: 5 -CGTC 
AACTTGCCGCCGTC-3' (modified from [24]); histone 
H4 [LinJ35.1400, GenBank: XM_001468907] forward: 5'- 
ACACCGAGTATGCG-3 ' and reverse: 5 -TAGCCGTA 
GAGGATG-3' [9]; SHERP [GenBank: XM_003392466] 
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forward: 5 -CAATGCGCACAACAAGATCCAG-3 ' and 
reverse: 5 -TACGAGCCGCCGCTTATCTTGTC-3 ' [9]; 
rRNA4S [GenBank: CC144545] forward: 5-CCTACCA 
TGCCGTGTCCTTCTA-3 ' and reverse: 5-AACGACCC 
CTGCAGCAATAC-3 ' [25]. Changes in relative gene 
expression were determined with AACT method and re- 
sults show fold changes comparative to day 1 calculated 
by2' AACT 

In vitro infections 

Bone marrow-derived macrophages (BMMo) were pro- 
duced as described previously [9]. Stationary promas- 
tigotes cultivated in NNN for 4 days were washed and put 
in contact with the cells in 1:10 ratio (celbparasites) for 
4 hours. Extracellular parasites were washed away with 
PBS and the cells incubated for more 24, 48, 72 or 
96 hours or fixed immediately with 2% PFA. The macro- 
phages were mounted on Vectashield with DAPI (Vector 
Laboratories) and 100 infected cells or 400 total cells were 
counted in duplicate by fluorescence microscopy in a 
Zeiss Axioskop (Carl Zeiss). The percentage of infected 
cells and the geometric mean of the number of parasites 
per infected cell were evaluated. The infection index was 
calculated by multiplication of both parameters to account 
for the overall parasite load. 

In vivo infections 

7-8 week-old BALB/c male mice (4-5 animals per 
group, except in 2-week infections with E390M strain 
where only 3 animals were used) were infected via the 
intraperitoneal route with 10 8 promastigotes of each strain 
cultivated in NNN for four days. After 2 or 6 weeks of 
infection mice were anesthetized with isoflurane and 
sacrificed by cervical dislocation. Blood, inguinal lymph 
nodes, spleen, liver and femoral bone marrow were recov- 
ered for quantification of parasite load. Blood and spleen 
were also used for the evaluation of humoral and cellular 
responses. 

Parasite load quantification 

Parasite load was quantified in samples that were col- 
lected and frozen at the time of animal sacrifice. We 
used 200 uL of blood collected with EDTA, 10 mg of 
spleen and liver (single cell suspensions), 3 x 10 6 bone 
marrow cells and the inguinal draining lymph node to 
extract DNA. First, 400 uL of a buffer containing 10 mM 
NaCl, 10 mM EDTA and 10 mM Tris-HCl with pH 8.0 
were added to the samples, which were incubated over- 
night with 40 ug of proteinase K ( Sigma- Aldrich) at 56°C 
with shaking. Then, the samples were vortexed and 
incubated for 20 minutes at 70°C. DNA was extracted 
using phenol/chloroform/isoamyl alcohol (all from Merck 
Millipore). After precipitation with ice-cold 70% ethanol 
solution, DNA was dissolved in 100-200 uL of nuclease- 



free water. We quantified the total DNA in a Nanodrop 
ND-1000 spectrophotometer and prepared dilutions of 
concentrations adjusted for each tissue. We quantified 
Leishmania sp. DNA by qPCR using 1000 nM of R223 
and 500 nM of R333 primers (Sigma- Aldrich) for the 
small subunit rRNA (SSUrRNA) [26]. Depending on the 
tissue, 100 to 400 ng of total DNA served as a template in 
a 20 uL reaction using LightCycler FastStart DNA Master 
SYBR Green I kit (Roche Applied Science) according to 
the manufacturers instructions, in a touchdown qPCR 
performed in a LightCycler 2.0 carousel-based instrument 
(Roche Applied Science) with final annealing temperature 
of 65°C [27]. CTs were extrapolated in a standard curve 
constructed with serial dilutions of L. infantum DNA 
(strain JPC, MCAN/ES/98/LLM-722) diluted in host 
DNA (from spleen of naive mice) to calculate Leishmania 
content in parasites/ ug DNA. Whenever the qPCR gave a 
positive (with the expected melting curve) but unquantifi- 
able value or a doubtable specific product (aberrant melt- 
ing curve), we performed a nested PCR [28] that has a 
higher sensitivity (0.01 parasites) than the qPCR (0.6 para- 
sites) to confirm the positivity of the quantitative result. 
300 nM of R221 and R332 primers [26] were used for the 
first amplification reaction. For the second reaction, 10 uL 
of the first PCR product diluted 1:40 served as template 
with the same R223 and R333 primers (300 nM and 150 
nM, respectively) used for the qPCR. This molecular 
quantification was applied after proper validation by com- 
parison with limiting dilution assay (Additional file 1: 
Additional Methods and Additional file 2: Figure SI). 

Splenic cell populations 

5 x 10 5 splenocytes were surface-stained for 20 minutes 
at 4°C with saturating concentrations of monoclonal 
antibodies (all from Biolegend). After washing twice with 
PBS/FBS 2%, the cells were examined by flow cytometry 
in a FACSCanto (BD Bioscences) and analyzed with 
Flowjo software. After acquisition of 50000 cells identi- 
fied by FSC and SSC parameters, major populations 
were identified as follows: CD4 + T lymphocytes (PerCp. 
Cy5.5 anti-CD3, clone 17A2; APC.Cy7 anti-CD4, clone 
GK 1.5), CD8 + T lymphocytes (PerCp.Cy5.5 anti-CD3; 
FITC anti-CD8, clone 53-6.7), B cells (FITC anti-CD19, 
clone 6D5), monocytes/macrophages (PE.Cy7 anti-CD lib, 
clone Ml/70; PerCp.Cy5.5 anti-Ly6C, clone HK1.4). 

Leishmania-spec\f\c immunoglobulins 

The specific humoral response was analyzed by ELISA 
as described elsewhere [18]. In short, 96-well microtitration 
plates (Greiner Bio-One) were coated with 10 ug/mL of 
soluble Leishmania antigens (SLA) in carbonates buffer 
pH 8.5 and then blocked with PBS/gelatin 1%. Sera were 
diluted 1:100 and incubated for 2 hours at 37°C. After 
washing with PBS/tween 20 0.1%, HRP-conjugated anti- 
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IgGl or anti-IgG2a (Southern Biotech) were added to the 
wells at a dilution of 1:5000 and incubated for 30 minutes 
at 37°C The plates were revealed with 0.5 mg/mL of o- 
phenylenediamine dihydrochloride (Sigma- Aldrich) in cit- 
rate buffer pH 4.0 and the reaction was stopped with HQ 
3 N. The absorbance was read at 492 nm in a Synergy 2 
microplate reader (Biotek). 

Animals and ethics statement 

For the in vitro experiments we used 10-12 week-old 
female BALB/c mice bred and maintained at IBMC - 
Instituto de Biologia Molecular e Celular (Portugal) ani- 
mal facilities. For the in vivo experiments 7-8 week-old 
male BALB/c mice were bred and maintained at the 
Instituto de Salud Carlos III (Spain) animal facilities. 
Mice were housed in IVC cabinets with sterile food and 
water ad libitum. All experiments conducted were car- 
ried out in accordance with the IBMC.INEB and ISCIII 
Animal Ethics Committees and the Portuguese and 
Spanish National Authorities for Animal Health guide- 
lines that follow the statements on the directive 2010/ 
63/EU of the European Parliament and of the Council. 
ACS has an accreditation for animal research given from 
Portuguese Veterinary Direction (Ministerial Directive 
1005/92). 

Statistical analysis 

GraphPad Prism 5 (GraphPad Software) was used to 
perform all the statistical analysis. The results are 
presented as means ± standard deviations (SD). To com- 
pare statistical differences between means two-sided t 
test or one-way ANOVA followed by Dunnett s multiple 
comparison test were run when comparing 2 or more 
groups, respectively, unless otherwise stated. * p < 0.05, 
** p < 0.01 and *** p < 0.001. 

Results and discussion 

Molecular characterization of the clinical isolates of L 
infantum 

To understand the intraspecies polymorphisms and its 
possible impact on both in vitro and in vivo infectivity, 
we characterized certain molecular aspects of these four 
L. infantum strains. 

A molecular approach was followed by us aiming not 
only to confirm the identity [19] of Bibiano, E390M and 
HL strains, together with the laboratory's standard L. 
infantum (ST) strain, but also to subtype the isolates 
according to Haralambous et al [21] and Kuhls et al [20]. 
Molecular genotyping of the strains indicated that all four 
were 100% consistent with L. infantum and were clustered 
in the same ITS type A group (data not shown), which is 
the most common in specimens from the Mediterranean 
area even within different zymodemes [20]. Moreover, 
Bibiano, HL and ST were classified in the I<26 group lb, 



the most frequent in the Iberian Peninsula [21], while 
E390M belongs to a new I<26 group reported here for the 
first time, since it returned an 836 bp amplicon [GenBank: 
KC576808] never found before (Figure 1A and Additional 
file 3: Figure S2). 

PCR-RFLP of kinetoplast DNA minicircles was used as 
a tool for creating an individual identity for each strain. 
Because we were working with L. infantum strains that 
preferably could have similar growth and morphology, 
they would be indistinguishable in in vitro cultures. In 
case of cross-contaminations between strains [29,30], we 
would like to have a tool for identification of our parasites. 
All the strains showed very complex profiles (Figure IB), 
but each one of the specimens was clearly identifiable by 
the examination of the most intense bands. 

Characterization of biological features of the 
promastigotes generated in vitro 

In vitro cultivation has a major impact on the virulence 
of the pathogens due to intrinsic properties of the culture 
media that modulate Leishmania infectivity (Santarem, N. 
and Cunha, J., submitted results and [8,31]), or to the loss 
of adaptive capacities to mammalian host cells resulting 
from long-term in vitro cultivation of promastigotes. A 
relevant experimental bias can be introduced if these fac- 
tors are not considered. Hence, the comparison of the in- 
fective capacity of distinct strains should take into account 
the adaptation to culture conditions and/or the axenic 
growth behavior [7,32] . 

We started to characterize the axenic growth of the 
four L. infantum strains. Aiming to identify any differ- 
ence in their infectivity that could justify the dissimilar 
outcomes of the disease in the natural infections, we first 
had to understand what specific nutritional needs the 
parasites might have in order to standardize the in vitro 
culture conditions. Preliminary experiments with para- 
sites cultivated in RPMI medium exposed remarkable 
differences among the growth of each L. infantum strain 
(data not shown). Attempts to find a culture medium 
suitable for the continuous and comparable growth of 
the four strains included doubling the FBS content to 
20% and adding glucose to RPMI medium, the use of 
Schneider s Insect medium and the preparation of a mix- 
ture of RPMI and Schneider media in equal parts (data 
not shown; the composition of these culture media is 
detailed in Additional file 1: Additional Methods). How- 
ever, only the use of Novy-MacNeal-Nicolle medium 
(NNN) allowed us to guarantee an in vitro homogeneous 
culture condition where all four strains could be easily 
cultivated, maintained and have a similar development 
that would allow the direct comparison of strain infect- 
ivity. Bibiano, E390M, HL and ST strains showed per- 
fectly overlapped growth curves in NNN (Figure 2A), 
with morphologically indistinguishable promastigotes 
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Figure 1 Molecular characterization of L infantum isolates. (A) Amplification of the k26 gene. See Additional file 3: Figure S2 for coding DNA 

sequences and alignment. (B) RFLP patterns for L infantum isolates after Hoe\\\ digestion of an amplified fragment of the minicircles kDNA. 

Profiles for DNA extracted from axenic promastigotes are shown, though DNA from experimentally infected murine tissues delivered comparable 

restriction patterns. For each strain, left lane corresponds to the digested DNA and right lane to undigested product. Gels were prepared with 2% 

(A) or 2.5% (B) agarose and stained with ethidium bromide. MW - molecular weight marker with bp units assigned in the left of the figure, 
k J 



that maintained high viability (90-95%) at least until the 
fourth day of culture (Figure 2B). 

The culture of all promastigote strains in NNN en- 
abled a rapid exponential parasite growth that peaked at 
day 3 (~10 8 parasites/mL) after which it stabilized until 



day 6. These timings could be confirmed by the analysis 
of the cell cycle (Figure 2C and Additional file 4: Figure 
S3). On the first day of culture, around 50% of the para- 
sites were in S/G2 phases in a clear indication of active 
proliferation. From the third day on, only « 20% of the 




1 0 4 "l 1 1 1 1 1 1 O-l 1 1 1 1 1 r 

0123456 0123456 
Days of culture Days of culture 



C D 




Days of culture 

Figure 2 In vitro development of L infantum strains. (A) Growth curves of Bibiano (light grey), E390M (dark grey), HL (black) and ST (white) L. 
infantum strains in NNN. Cultures were started with 10 6 promastigotes of each strain per ml_ and growth was followed for a week, with a daily 
counting of parasite numbers in a Neubauer chamber. (B) Parasite viability measured by flow cytometry as percentage of AnnV77AAD" cells. The 
mean of three independent experiments is plotted; bars represent SD. (C) Cell cycle analysis of HL during growth period. Data show means of 
one representative experiment of three independent assays. (D) Illustrative phenotype of L infantum promastigotes with 4 days of culture in NNN 
captured by flow cytometry. 
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cells were found to be in division. Moreover, during the 
stationary phase promastigotes with small bodies and 
large flagella could be clearly identified by optic micros- 
copy (data not shown); these were confirmed by flow cy- 
tometry (Figure 2D) as FSC lo Pr cells [33,34]. Based on 
this typical morphology of metacyclic promastigotes, we 
analyzed metacyclogenesis throughout the 6 days of cul- 
ture by quantifying the expression of specific genes which 
are upregulated (Meta-1 and SHERP) or downregulated 
(histone H4) during the process [25,33,35] (Figure 3 and 
Additional file 5: Figure S4). Despite the observed inter- 
strain variations on the fold modifications of each gene, all 
of the strains showed a dramatic decrease in histone H4 
expression consistent with the cell cycle analysis. For 
Bibiano, Meta-1 gene expression significantly increased 



every day of culture, while for E390M the major fold 
changes were detected in the SHERP gene expression. 
Both Meta-1 and SHERP increased overtime for HL and 
ST until day 4 and after that they recovered to levels 
close to those of day 1. We cannot directly compare the 
metacyclogenesis process between the four strains, but 
we must take into account the total information avail- 
able to affirm that these promastigotes were differentiat- 
ing into metacyclic forms. The leading factor is that the 
overall analysis of the three genes studied pointed to- 
wards metacyclogenesis. This trend could be better 
understood when calculating the ratio between the up 
and downregulated genes (Additional file 5: Figure S4), 
as all the strains presented an increase in the ratios over 
the 6 days of culture. 
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Figure 3 Indirect measurement of metacyclogenesis. Quantification of (A) Meta-1, (B) SHERP and (C) histone H4 transcription by RT-PCR over 
the time of culture in NNN medium for the four strains. Bars represent the mean fold change relative to day 1 with SD of two independent 
experiments. Statistical significant differences between day 1 and the following days were determined with One-way ANOVA and Dunnett's 
multiple comparison test. 
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Evaluation of the infectivity of L infantum isolates 

To evaluate parasite infectivity, we infected macrophages 
and followed the progression of the infection for four 
days (Figure 4). At the initial time-point, the four strains 
infected approximately all the cells without any differ- 
ences in the total percentage of infected macrophages 
(Figure 4A). However, when evaluating the parasite load 
in each infected cell, HL promastigotes were more ef- 
fective in the invasion of the macrophages with a mean 
of 12 ± 1.5 parasites in each infected cell versus 3 to 5 par- 
asites found for the other strains (Figure 4B). After the 
first 24 hours, an accentuated reduction in the infection 



index was detected for the four strains (Figure 4C), which 
was maintained, though softened, during the entire period 
of the assay. Nevertheless, HL parasites were more resist- 
ant to macrophage-specific killing machinery as displayed 
by the increased infection index throughout the study. For 
this reason, this strain was considered to be more infective 
than the ST virulent strain [36,37]. Interestingly, 72 and 
96 hours post-infection, macrophages could harbor sig- 
nificantly more Bibiano than ST parasites, though it did 
not translate into a higher overall parasite load as shown 
for HL (Figure 4C) because of the low percentage of 
infected cells. 
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As is well known, disease manifestation depends on 
the virulence of Leishmania strain [12,17] but also on 
the genetic background [38] and immune status of the 
host [14,39]. Hence, the differences in infectivity detected 
in vitro were explored in the BALB/c model of visceral 
leishmaniasis to evaluate the influence of the parasites' in- 
trinsic characteristics in the ability to cause the disease. 
We therefore studied parasitological and immunological 
features in the acute and chronic phases of infection. Para- 
site load was quantified in the spleen, liver and bone mar- 
row as the main target organs in this model. The presence 
of Leishmania in the draining lymph nodes and the blood 
were also investigated (Figure 5). 

In the acute phase, 2 weeks after the infection, no 
significant differences were calculated between the four 
strains in all the five tissues studied. Despite the lack of 
statistical significance, important differences could be 
appreciated between strains. In the visceral organs 
(Figures 5A and B), HL and Bibiano showed very similar 
parasite burden, representing the most infective strains 
in these tissues. As compared to ST, they presented over 



1 logarithm higher parasite load in spleen and were ~ 
800 times higher in the liver. ST was indeed the least effi- 
cient strain colonizing these organs, as E390M was able to 
infect almost 3 and 15 times more in the spleen and liver 
of mice. In the inguinal lymph nodes (Figure 5D) and 
the blood (Figure 5E), HL was the strain that showed 
higher parasitism. Its values surmounted 14, 300 and 
1100 times the burden of Bibiano, E390M and ST, re- 
spectively, in the lymph nodes, and 3, 20 and 46 times, 
respectively, in the blood. The bone marrow was more 
infected by ST, though without accentuated differences 
between the other strains (Figure 5C). After this ana- 
lysis one could expect ST to be the least infective strain 
in the acute phase of the disease, once HL was the 
strain that showed higher infectivity in the majority of 
the tissues examined, followed close by Bibiano. None- 
theless, estimating the total number of parasites existing 
in the whole animal, we calculated « 5 x 10 8 ST para- 
sites, approximately double the number estimated for 
Bibiano and HL and over a logarithm more than E390M 
(Table 1). 
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Figure 5 Quantitative distribution of L. infantum in BALB/c mice 2 and 6 weeks after infection. The parasite burden was quantified on the 
(A) spleen, (B) liver, (C) bone marrow, (D) lymph node and (E) blood. Data represent means ± SD of 3-5 animals per group of a representative 
experiment of two independent assays. T test was used to determine statistically significant differences between 2 and 6 weeks of infection with 
each strain; * p < 0.05. With one-way ANOVA followed by Tukey's multiple comparison test we calculated statistically significant differences 
between strains at each time point; # p < 0.05. Dashed line indicate limit of detection for quantification for each tissue. 
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In the chronic phase, 6 weeks after infection, the four 
strains showed relative infection profiles in the spleen 
(Figure 5A), bone marrow (Figure 5C) and lymph nodes 
(Figure 5D) similar to the acute phase. As before, HL 
was shown to be the most infective strain in spleen and 
lymph nodes. In the bone marrow the differences were 
once more not as accentuated between strains, although 
the higher parasite loads were found in HL infected 
mice. In the liver (Figure 5B), HL was still the most in- 
fective strain, but Bibiano, which in the other tissues 
presented comparable parasite loads, was, significantly, 
the least infective strain, whereas E390M and ST pro- 
duced intermediate infections. In the blood (Figure 5E), 
the differences observed at 2 weeks post-infection were 
neutralized, as the four strains showed similar levels of 
circulating parasites. 

Distribution and compartmentalization throughout the 
infection 

Evaluating the progression of the disease, these four 
strains depicted very different trends. In the acute phase 
of infection, Bibiano and HL were found in very high 
numbers in the visceral organs but evolved in the oppos- 
ite directions with chronicity. Bibiano was efficiently 
cleared from the liver, with a 425-fold reduction, though 
in the spleen the parasite load did not alter. As to HL, 
the liver infection was maintained and the splenic para- 
site burden increased 5 times from 2 to 6 weeks after in- 
fection, which showed not only high capacity to infect 
but also to perpetuate in the host. E390M, on the con- 
trary, showed a low infective phenotype, with the lowest 
parasite loads in all the tissues quantified in the acute 
phase of infection. Through time, this strain was not 
able to proliferate in the spleen or in the bone marrow; 
the parasites resisted in the liver and showed a 7.4-fold 
increase in the blood. Despite with disparate initial para- 
site loads (more than 6-fold difference), Bibiano and 
E390M, both isolated from HIV+ patients, followed a 
very similar trajectory in the progression of the disease 



and, eventually, would be eliminated over time in these 
immunocompetent BALB/c mice. The high parasitemia 
presented by these two strains in the chronic phase may 
facilitate the anthroponotic transmission of Leishmania 
between the intravenous drug users (IVDUs) [42], one of 
the populations with highest risk of WW I Leishmania co- 
infection [43]. Concerning ST, the standard virulent 
strain in our laboratory, showed a clear tropism for the 
bone marrow in the acute phase, with the lowest parasite 
loads in the remaining organs compared to other strains. 
However, 6 weeks post-infection, ST dramatically multi- 
plied reaching levels ~ 6- and ~ 17-fold higher in the 
spleen and liver, respectively. 

As a final remark, this study allowed us to verify that 
bone marrow parasite load is maintained in a range that 
does not suffer major alterations either over- time nor is 
it strain-dependent. Possibly this is the reason why bone 
marrow aspirates are the eligible sample for leishmania- 
sis diagnosis, either for microscopic analysis, culture or 
molecular techniques [4]. 

Infectivity relates to cell modulation 

Along with the differences described in the parasite load, 
we analyzed the weight of spleen and liver of the 
infected animals and compared them with naive mice 
(Additional file 6: Figure S5A and B). In general, small 
fluctuations were detected in the weight of the organs, 
which followed the trends exposed above for the parasite 
loads. Indeed, HL infected mice presented hepatos- 
plenomegaly in the chronic phase of the disease, the 
pathologic hallmark of visceral leishmaniasis (reviewed 
in [44]). We attributed the splenomegaly not only to the 
presence of elevated numbers of Leishmania but also to 
the expansion of the main cellular populations detected 
in the spleen (Figure 6). A significant increase in both 
CD4 + and CD8 + T cells, B cells and macrophages was 
quantified in mice infected for 6 weeks with HL, an incre- 
ment that corresponded to double the numbers found in 
the naive mice. While the expansion of lymphocyte 



Table 1 Estimated overall parasite load of L. infantum-'infected BALB/c mice 2 and 6 weeks post-infection 

Total parasite load 2 weeks post-infection 6 weeks post-infection 



Bibiano E390M HL ST Bibiano E390M HL ST 



Spleen 


67.7 x 10 4 


8.82 X10 4 


44.1 x 10 4 


3.81 x 10 4 


62.7 X10 4 


0.40 X10 4 


329 X10 4 


18.9x 10 4 


Liver 


157 x 10 6 


2.77 X10 6 


1 1 2 x 1 0 6 


0.14X 10 6 


0.35 X10 6 


2.54 X10 6 


86.8 X10 6 


2.42 X10 6 


Bone marrow 


15.0X 10 7 


4.49 X10 7 


11.9x 10 7 


50.9 x 10 7 


13.9X10 7 


2.29 X10 7 


18.1 X10 7 


11.7X 10 7 


Lymph nodes 


4.83 x 10 4 


0.19X10 4 


67.0 x 10 4 


0.1 Ox 10 4 


1.78X10 4 


0.015 X10 4 


10.3 x10 4 


0.59 X10 4 


Blood 


0.97 x 10 3 


0.42 X10 3 


19.4X 10 3 


0.32 X10 3 


6.54 X10 3 


3.12X10 3 


1 .95 x 1 0 3 


4.53 x 10 3 


Whole animal 


3.08 x10 s 


0.48 x10 s 


2.32 x10 s 


5.09 x 10 s 


1.40x10 s 


0.25 x10 s 


2.71 X10 8 


1.19X 10 s 



Mean weights of spleen, liver and inguinal lymph nodes of the infected animals (Additional file 6: Figure S5) were used in these calculations. Values for total 
volume of blood (1.8 mL) and total number of bone marrow cells (4.5 x 10 8 ) were estimated as published elsewhere [40,41]. 
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Figure 6 Cell populations in spleens of naive and Le/sftman/a-infected mice in the acute and chronic phases. (A), (F) Total cells were 
counted and stained for identification of major splenic populations by flow cytometry. (B), (G) CD4 + T cells. (C), (H) CD8 + T cells. (D), (I) B cells. 
(E), (J) Monocytes/macrophages. (A)-(E) Absolute number of cells per spleen. Data represent means ± SD of 3 to 5 animals per group of one 
experiment representative of two. (F)-(J) Fold modification of cell numbers in the infected mice in relation to naive. Boxes and whiskers with 
5-95 percentile and mean (showed by "+") of 3 to 5 animals. Dashed and solid lines indicate a 2- or 4-fold modification, respectively, relative to 
naive mice. One-way ANOVA followed by Dunnett's multiple comparison test was performed to calculate statistically significant differences 
between naive and infected mice at 2 and 6 weeks after infection. 



populations are classically linked to cellular and humoral 
immune responses, the presence and increase of CDllb 
+ Ly6C + monocytes were recently described to play a 
major role in the architectural remodeling of the spleen 
during experimental visceral leishmaniasis, mainly in 
the vascularization of the red pulp that accompanies 
splenomegaly [45]. 



The enlargement in the B cell population explains the 
elevated titers of anti- Leishmania antibodies quantified 
in the chronic infection by HL (Figure 7). Both IgG2a 
and IgGl were generated in high levels, leading to the 
exacerbation of the pathology, as described by others 
[46-48]. However, in E390M infections we did not detect 
any B cell expansion, as measured at 2 weeks post 
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Figure 7 Leishmania-spedflc humoral response. Leishmania-speclflc sera reactivity of naive and infected animals 2 and 6 weeks post-infection 
was analyzed by ELISA. Specific (A) lgG2a and (B) lgG1 were quantified and are depicted as means ±SD of one representative experiment of 2 
independents. Statistically significant differences are pointed out as given by one-way ANOVA followed by Dunnett's multiple comparison test. 



infection and later at 6 weeks, though IgG2a and IgGl 
were significantly increased compared to age-matched 
naive mice. We speculate that this antibody production 
could be in part related to the I<26 gene. Leishmania k26 
(also known as HASPB) protein, as well as SHERP that 
share the same locus on chromosome 23 [49], are stage- 
regulated proteins, expressed only in the mammalian host 
infective forms (metacyclic promastigotes and amastigotes) 
[50]. k26 has a central core composed of a 10-11 amino 
acid repeats (PKEDGHTQKND/PKEDGRTQKN in Add- 
itional file 3: Figure S2) which present high inter and 
intra-specific variability [49,51,52], a feature that makes 
it an interesting tool for molecular typing of different 
Leishmania species and strains [21]. Other than this, 
k26 is a highly immunogenic antigen [49] with proven 
efficacy as a vaccine in murine models of visceral leish- 
maniasis by L. donovani [51,53] and partial efficacy in 
canine leishmaniasis by L. infantum [54]. Its immuno- 
genic properties make k26 an interesting antigen that 
has been studied for diagnostic purposes [55-57]. As we 
reported in this work, E390M k26 protein has more re- 
peats than MON-1 strains, which may influence the 
type and the strength of the humoral response as those 
amino acid repeats were determined to be B cell epi- 
topes [52]. We point towards this argument since 
Bibiano and ST strains, that proved to be more infective 
than E390M, were not able to produce specific anti- 
bodies nor increased splenic cellularity in the acute or 
the chronic phases of murine VL. 

Conclusions 

The molecular typing strategy confirmed the previous 
zymodeme characterization by MLEE and provided fur- 
ther knowledge that can be applied to diagnostics and 
population genetics studies. In this sense, the I<26 se- 
quence for E390M strain generated in this work adds in- 
formation for a zymodeme (MON-284) not included in 
the study by Harambolous [21], thus contributing to k26 
gene-based typing methodology, which is being used 



increasingly in population genetics and molecular epi- 
demiology studies related to the L. donovani complex 
[58,59]. Laboratory conditions for in vitro culture were 
set to produce Bibiano, E390M, HL and ST fit pro- 
mastigotes in the same developmental stage. In vivo in- 
fections with HL confirmed the in vitro phenotype of 
the most infective strain as more parasites were esti- 
mated to be present in the whole animal. ST was also 
considered to be highly infective though with a slower 
progression over time. Bibiano and E390M, isolated 
from HIV+ patients, showed differential infectivity and 
immunomodulation that could be influenced by the ini- 
tial compartmentalization in host tissues. Interestingly, 
the most (HL) and the least (E390M) infective strains were 
the most immunogenic, revealing high levels of anti-Leish- 
mania IgG2a and IgGl, especially in the chronic phase of 
infection. 

This work is in line with previous studies [6,7,17,60] 
that show that leishmaniasis is a multifactorial disease 
and the broad spectrum of clinical manifestations de- 
pends on the genetics and inherent characteristics of 
the parasite coordinated with the susceptibility of the 
host. 



Additional files 



Additional file 1: Additional methods. 

Additional file 2: Figure S2. Validation of the qPCR methodology for 
quantification of the parasite loads in murine tissues. BALB/c mice were 
infected for 2 weeks with 1 0 s /_. infantum strains (3 mice for each strain) 
cultivated in an equivolumetric mixture of RPMI and Schneider for 5 days. 
(A) Spleen and (B) liver were collected and cell suspensions were 
prepared. Parasite loads were quantified by qPCR and compared with the 
values obtained by limiting dilution assay. 

Additional file 3: Figure S2. k26 gene alignment. Coding DNA 
sequences (cds) were translated into amino acids and aligned with 
standard sequences generated by Haralambous et al. [21] [GenBank: 
EF504255 and EF504256]. Complete cds of k26 gene [GenBank: 
XM_00 1465758.2] is aligned in the first row. The four L. infantum isolates 
are aligned in the bottom rows. E390M k26 sequence was entered in the 
GenBank with the access number KC576808. 
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Additional file 4: Figure S3. Cell cycle analysis. (A) Representative 
histograms of PI content at day 1 and 6 of promastigote culture. (B-D) 
The percentage of (B) Bibiano, (C) E390M and (D) ST promastigotes in G1, 
S or G2 phases of the cell cycle was determined at days 1 to 6 of culture 
in NNN medium by flow cytometry after staining the parasites with a PI 
solution. Three independent experiments were performed, one 
representative experiment is shown. 

Additional file 5: Figure S4. Variation of the transcription of 
metacyclogenesis-dependent genes. Meta-1, 5HERP and histone H4 
transcription was quantified by RT-PCR over the time of culture in NNN 
medium. To evaluate the progression of metacyclogenesis we related the 
up-regulation of (A) Meta-1 and (B) SHERP with the down-regulation of 
histone H4 applying a mathematical ratio. Bars represent the mean fold 
change relative to day 1 with SD of two independent experiments. 
Statistically significant differences between day 1 and the following days 
were determined with One-way ANOVA and Dunnett's multiple 
comparison test. 

Additional file 6: Figure S5. Organ weight 2 and 6 weeks post- 
infection. Mice were sacrificed and (A) spleen, (B) liver and (C) inguinal 
lymph nodes were collected and weighted. Data represent means ±SD 
of 3-5 animals of one experiment representative of two. One-way 
ANOVA followed by Dunnett's multiple comparison test were run for 
statistical analysis between groups in each time point. 
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